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Abstract 

Background: Lipase from Rhizopus chinensis is a versatile biocatalyst for various bioconversions and has been 
expressed at high-level in Pichia pastoris. However, the use of R. chinensis lipase in industrial applications is restricted 
by its low thermostability. Directed evolution has been proven to be a powerful and efficient protein engineering 
tool for improvement of biocatalysts. The present work describes improvement of the thermostability of R. chinensis 
lipase by directed evolution using P. pastoris as the host. 

Results: An efficient, fast and highly simplified method was developed to create a mutant gene library in P. pastoris 
based on in vivo recombination, whose recombination efficiency could reach 2.3 x 10 s /ug DNA. The thermostability 
of r27RCL was improved significantly by two rounds of error-prone PCR and two rounds of DNA shuffling in 
P. pastoris. The S4-3 variant was found to be the most thermostable lipase, under the conditions tested. Compared 
with the parent, the optimum temperature of S4-3 was two degrees higher, T" m was 22 degrees higher and 
half-lives at 60°C and 65°C were 46- and 23- times longer. Moreover, the catalytic efficiency k cat /K m of S4-3 was 
comparable to the parent. Stabilizing mutations probably increased thermostability by increasing the hydrophilicity 
and polarity of the protein surface and creating hydrophobic contacts inside the protein. 

Conclusions: P. pastoris was shown to be a valuable cell factory to improve thermostability of enzymes by directed 
evolution and it also could be used for improving other properties of enzymes. In this study, by using P. pastoris as 
a host to build mutant pool, we succeeded in obtaining a thermostable variant S4-3 without compromising 
enzyme activity and making it a highly promising candidate for future applications at high temperatures. 
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Background 

Lipases (triacylglycerol ester hydrolases EC 3.1.1.3) are 
enzymes that hydrolyze the ester bonds of water- 
insoluble substrates at the interface between substrate 
and water. They are well known because of their re- 
markable levels of activity and stability in non-aqueous 
environments, which makes them as attractive enzyme 
for use in industrial applications, such as in food, 
pharmaceutical, chemical and detergent industries [1-3]. 

Rhizopus lipases have been widely used in food indus- 
tries [4,5]. Nevertheless, most Rhizopus lipases are meso- 
zymes as they are produced by mesophilic organisms 

* Correspondence: yxu@jiangnan.edu.cn 
+ Equal contributors 

'State Key Laboratory of Food Science and Technology, Jiangnan University, 
1800 Lihu Avenue, Wuxi 214122, China 

2 Key Laboratory of Industrial Biotechnology, Ministry of Education, School of 
Biotechnology, Jiangnan University, 1800 Lihu Avenue, Wuxi 214122, China 
Full list of author information is available at the end of the article 

Bio Med Central 



and exhibit unfavorable thermostability [6]. For example, 
lipase RAL from R. arrhizus showed poor thermostabil- 
ity when incubated at 30°C for 30 min [7]. Lipase ROL 
from R. oryzae showed highest activity at 35°C and about 
35% of the lipase activity was lost after incubation at 45° 
C for 30 min [8]. Another lipase proROL derived from 
R. oryzae DSM853 exhibited better thermostability with 
a T m value of 57.5°C [9]. In our previous study, a lipase 
gene proRCL from R. chinensis CCTCC M201021 was 
cloned and expressed at high-level in Pichia pastoris as 
r27RCL, whose activity was about 580 times higher than 
that of wild-type R. chinensis lipase [10]. This recombin- 
ant enzyme showed a high potential for industrial usage, 
as an additive for bread baking, synthesis of eicosapenta- 
enoic acid (EPA), docosahexaenoic acid (DHA), sorbitan 
oleate, and ethyl esters [11-15]. However, the use of lip- 
ase r27RCL from R. chinensis is restricted by its low 
thermostability. 
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Directed evolution has been proven to be a powerful 
and efficient protein engineering tool for improvement 
of biocatalysts [16]. Several lipases from Rhizopus sp. 
were subjected to directed evolution in an effort to pro- 
duce thermostable variants [7,17]. The half-life of RAL 
expressed in P. pastoris was prolonged by 12-fold at 
50°C [7]. The T m value of RNL from R. niveus 
expressed in Saccharomyces cerevisiae was improved 
from 47°C to 57°C [17]. 

The production of an active Rhizopus lipase has been 
performed in Escherichia coli [18], S. cerevisiae [19] and P. 
pastoris [20-22], among which the expression of Rhizopus 
lipase in P. pastoris was the best. The P. pastoris expres- 
sion system offers several advantages, as a desirable host 
in directed evolution, for easy library screening and fur- 
ther industrial application, such as the ability to perform 
eukaryotic protein modifications, high levels of protein ex- 
pression at the intra- or extracellular level and an alcohol 
oxidase 1 {AOX1) gene promoter tightly regulated by 
methanol [23,24]. However, mutant library construc- 
tion in P. pastoris usually requires time-consuming 
multistep procedures which include an intermediate 
bacterial host {E. coli) and restriction digestion pro- 
cedures [25]. 

In this work, we developed a faster and simplified mu- 
tant library production method, based on which, the 
thermostability of the lipase r27RCL from R. chinensis 
CCTCC M201021 was improved by directed evolution. 

Results and discussion 

Library construction 

A novel approach for lipase mutant library construction 
in P. pastoris was established based on the formation of 
a recombination cassette in vivo (Figure 1). Three pairs 
of overlapping primers were used to generate the mutant 
gene (fragment B/C) and two different vector fragments 
(A and D). Primer BC-F had 66 bp overlapped pairs with 
the 5'-end of the gene sequence ("a" in Figure 1) of frag- 
ment D and primer BC-R had 60 bp overlapped pairs 
with the 3'-end of the gene sequence ("b" in Figure 1) of 
fragment A. Therefore, fragment B/C generated by pri- 
mers BC-F and BC-R had 66 bp and 60 bp overlapped 
sequences with vector fragments D and A. The lipase 
gene fragment B (or C) and vector fragments A and D 
were mixed and transformed into P. pastoris competent 
cells by electroporation. In vivo, the fragments A, B/C 
and D self-reassembled to form a recombination cassette 
relying on the 60 bp and 66 bp homologous arms be- 
tween them (Steps 1 and 2 in Figure 1). The self- 
reassembled expression cassette was integrated into the 
genome of P. pastoris by homologous recombination of 
the flanking AOX1 gene sequences (Steps 3 and 4 in 
Figure 1). 



Recombination efficiency is one of the most important 
factors in directed evolution. In order to improve the re- 
combination efficiency, the effects of the fragment molar 
ratio and the length of the overlapping sequences for re- 
combination were investigated. As shown in Table 1, the 
recombination efficiency was improved with the 
increased molar ratio of fragment B/C:A:D in the DNA 
mixture by electroporation. About 2.3 x 10 5 / ug DNA 
recombinants were obtained when the molar ratio 
reached 10:1:1 and the values dropped to 1.8x10 /ug 
DNA, 1.2 x 10 s /ug DNA or 1.0 x 10 s /ug DNA when the 
molar ratio was 5:1:1, 3:1:1 or 1:1:1. In order to obtain 
high recombination efficiency, the fragment molar ratio 
B/C:A:D of 10:1:1 was chosen. The recombination effi- 
ciency was low with shorter overlapping sequences 
(<55 bp) due to the difficulty in formation of a recom- 
bination cassette in vivo, while too long overlapping 
sequences (>70 bp) resulted in PCR technical difficulties 
with blurred and weak target bands when amplified with 
such long primers. Medium length of overlapping 
sequences (55-70 bp) generated satisfied recombination 
efficiency without PCR technical difficulties (Table 1). It 
is well known that homologous recombination is a com- 
mon characteristic of yeast. Baudin et al. [26] demon- 
strated that the amount of homology necessary to 
promote efficient recombination-mediated gene disrup- 
tion in S. cerevisiae is in the order of 30-50 bp. Olden- 
burg et al. [27] reported that 20 bp of homology is 
sufficient to generate the PCR-directed recombination 
in S. cerevisiae, at least under certain conditions (vector: 
insert 1:6). 

In this work, in order to make directed evolution prac- 
tical in the host organism P. pastoris an efficient, fast 
and highly simplified method was developed to create a 
gene mutant library in P. pastoris. This method went a 
crucial step beyond previous approaches [25,28-30] be- 
cause it completely discarded time-consuming multi- 
step procedures including construction of recombinant 
plasmid, transformation of an intermediate bacterial 
host (E. coli), plasmid isolation and linearization. In our 
approach, only three PCR amplification processes and 
an electroporation step were required. Fernandez et al. 
[31] also reported a method to produce a mutant library 
in P. pastoris, which reassembled a recombination cas- 
sette in vitro and produced about 6000 colonies per 
1 ug of linear DNA. The method developed in our re- 
search was more convenient and efficient since self- 
reassembly of the linear recombinant cassette occurred 
in vivo. 

Thermostable variants 

The thermostability of r27RCL was improved by com- 
bining two rounds of error-prone PCR and two rounds 
of DNA shuffling (Figure 2). Error-prone PCR was used 
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Figure 1 (See legend on next page.) 
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(See figure on previous page.) 

Figure 1 The novel approach of construction of the lipase mutant library in P. pastoris based on the formation of a recombination 
cassette in vivo. The hatched areas (a) and (b) were overlapped sequences between amplified fragments 



to introduce an average of 1 to 2 amino acid substitu- 
tions per proRCL gene and DNA shuffling could com- 
bine beneficial mutations as well as create new 
mutations. After creating the first-generation mutant li- 
brary using error-prone PCR, approximately 5000 clones 
were screened for increased thermostability at 45°C. 
Two variants (epl-1 and epl-26) were identified and 
confirmed to have longer half-lives at 60°C (Table 2). 
epl-1 (P168L, V329A) and epl-26 (S234F) with a total 
of three substitutions were chosen as templates for gen- 
erating the second-generation mutant library using 
error-prone PCR. Only 0.1% of the second generation 
mutant pool (5200 clones) screened at 50°C exhibited 
better performance than the first-generation parents. 
The identified variants ep2-3, ep2-4, ep2-8 and ep2-12 
were selected for DNA shuffling. It was interesting to 
note that ep2-8 (S234F, P168L, K219D) and ep2-12 
(S234F, V329A, K161R) both acquired mutations that 
were derived from epl-1 and epl-26 via in vivo DNA 
shuffling. From each generation of DNA shuffling, ap- 
proximately 10000 clones were screened for increased 
thermostability. The third-generation evolution pro- 
duced three variants (S3-1, S3-4, S3-18) whose thermo- 
stability slightly increased, whereas the last generation 
of evolution yielded three variants (S4-3, S4-13, S4-20) 
with significantly improved thermostability (Figure 2). 
S4-3 was the most thermostable variant. It contained 
five mutations (A129S, A230T, S234F, L180H, H317P) 
observed in previous generations plus one new muta- 
tion (T218S) generated by the fourth round of DNA 
shuffling. 

Thermostabilities of the evolved lipases 

Table 2 lists the ty% values of the parent and the evolved 
thermostable variants during four generations of direc- 
ted evolution. At 60°C the t V2 of the parent was 4 min. 
The variants from the first and the second generations 
had half-lives of thermal inactivation two- to five-fold 
longer than the parent. The half-lives of the three var- 
iants discovered from the third round of directed evolu- 
tion (DNA shuffling) were 10 to 14 times longer than 
that of the parent at 60°C. The last round of mutagenesis 
resulted in three variants that had the longest half-lives 
of thermal inactivation, 40 to 46 times longer than the 
parent at 60°C, among which, S4-3 was the most 
thermostable enzyme. The enzymatic properties of the 
purified S4-3 were then studied and its thermo- 
mechanism was explored. 



Initial activity vs thermostability of variants 

In order to determine the relationship between thermo- 
stability and enzyme activity, the initial activities of all 12 
variants were measured. As shown in Figure 3, five var- 
iants exhibited improved T m as well as higher activity. 
The activities of the other 7 variants were not affected 
much by increased thermostability, among which S3-1, 
with the lowest activity, still retained about 90% activity 
compared to the parent enzyme. These results provided 
strong evidence that improvement of thermostability and 
maintenance of enzyme activity could be simultaneously 
achieved, and improved thermostability did not occur at 
the expense of activity at low temperatures [32-34]. 

Enzymatic properties of the thermostable variant S4-3 

The most thermostable variant in the last generation, 
S4-3, was purified and its enzymatic properties were 
investigated. As shown in Figure 4, the residual activities 
of S4-3 and the parent were determined after incubation 
at different temperatures for 30 min. S4-3 displayed 
greatly- increased thermostability in the range of 20°C 
-60°C and retained 50% residual activity at 67°C, whereas 
the parent started to denature at 40°C and was com- 
pletely inactivated at 60°C after 30 min heat-treatment. 

The thermal inactivation of S4-3 and the parent were 
characterized by incubation at 60°C and 65°C (Figure 5). 
The thermal inactivation curves of the parent at 60°C 
and 65°C were much steeper than those of S4-3, and the 

Table 1 The effects of the fragment molar ratio of B/C:A: 
D and the length of the overlapping sequences on 
recombination efficiency 

Method Recombination 

efficiency c (x 1 0 5 , 
recombinants /ug DNA) 

Fragment B/C:A:D 10:1:1 2.30 

(molar ratio) a ~ 

3:1:1 1.20 

1:1:1 1.05 

Length of overlapping 10 0.10 

sequences (bp) b ~ ~ 

40 1.25 

55 1.75 

70 2.32 

85 2.35 

a The length of overlapping sequence was kept at 65 bp. 

b The fragment molar ratio of B/C:A:D was kept at 10:1:1. 

c Recombination efficiency is the number of recombinant cells (recombinants) 

generated by 1 ug of DNA fragments in a transformation procedure. 
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Figure 2 (See legend on next page.) 
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(See figure on previous page.) 

Figure 2 Lineage of variants and amino acid substitutions by directed evolution. Newly- introduced mutations in each generation are 
marked with asterisks 



half-lives of S4-3 at 60°C and 65°C were 23- and 46-fold 
higher than those of the parent. In addition, the 
optimum temperature of S4-3 was also increased by 2°C 
compared with that of the parent (data not shown). 

Kinetic parameters of the parent and S4-3 towards 
j?NPP were calculated and listed in Table 3. S4-3 exhib- 
ited a 32% decrease in K m whereas a 33% increase in k C!it , 
resulted in almost identical catalytic efficiency k cat /K m 
compared to the parent. 

Thermo-mechanism of S4-3 

The three-dimensional homology model of r27RCL 
shows that the five mutations A129S, A230T, S234F, 
L180H, T218S were all present on the surface of the en- 
zyme and far from the catalytic domain (S265-D324- 
H377) except for the mutation H317P (Figure 6). It is 
consistent with the notion that the initial steps in pro- 
tein unfolding during the irreversible thermal denatur- 
ation primarily involve surface-located parts of the 
protein [35,36]. 

The contributions of mutations A230T and A129S to- 
ward thermostability could be interpreted based on the 
structure model which located them at the end of a-helix 
and in a loop on the protein surface. The substitution of 
A230 for T230 and A129 for S129 may have increased the 
hydrophilicity and polarity of the protein, which enhanced 
the thermostability of the enzyme [37]. Mutation L180H 
located in a loop on the protein surface introduced a 
positively-charged residue which may increase protein sta- 
bility by creating a more favorable surface charge 

Table 2 Thermostability of enzyme variants from each 
generation 



Generation 


Variant 


t^/2 (min, at 60°C) 


Fold Improvement 


0 


parent 


4.0 


1.0 


I 


ep 1-1 


9.8 


2.4 




ep 1 -26 


8.7 


2.2 


2 


ep 2-3 


11.3 


2.8 




ep 2-4 


21.7 


5.4 




ep 2-8 


14.4 


3.6 




ep 2-12 


13.5 


3.4 


3 


S 3-1 


38.8 


9.7 




S 3-4 


58.9 


13.9 




S 3-18 


48.0 


12.0 


4 


S4-3 


184.0 


46.0 




S4-13 


168.3 


42.1 




S4-20 


159.5 


40.0 



distribution [38,39], as found in engineered thermostable 
cold shock proteins [40]. The mutation S234F was 
located in an a-helix on the protein surface. It was sur- 
prising that replacement of an uncharged hydrophilic 
residue for a nonpolar hydrophobic residue (S234F) at 
the protein surface increased protein stability but this 
effect has been observed previously [40,41]. The H317P 
mutation is buried inside the protein in a p-strand. 
There are two assumptions that could be made to ex- 
plain the improved thermostability. The first assumption 
was the addition of a hydrophobic residue. In the stabi- 
lized variant S4-3, hydrophilic histidine was replaced 
by hydrophobic proline, whose side-chain engaged in 
hydrophobic contacts with the side chain of Y293. This 
interaction was further improved by surrounding 
hydrophobic amino acids such as L137 and V341 and 
therefore enhanced protein stability. The alternative ex- 
planation was a reduction of the entropy of unfolding. 
The proline residue, with its pyrrolidine ring, had the 
lowest conformational entropy, thus the mutation histi- 
dine to proline should decrease the entropy of protein 
in the unfolded state and stabilize the protein [37]. 

It must be noted that the above-assumed thermo- 
mechanisms of S4-3 have been observed in many 
evolved thermostable enzymes and naturally-occuring 
thermophilic and hyperthermophilic proteins [35,42]. 
Determination of the crystal structure of S4-3 should 
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Figure 3 Initial activity vs T m of the parent and variants. The 

cross point of dashed lines and solid lines represent parent enzyme 
r27RCL. Variants exhibited higher activity and improved T m were 
squared with solid lines. The relative activity of the parent taking 
pNPP as substrate was set as 100% 
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Temperature ( °C) 

Figure 4 The effect of temperature on the stability of the 
parent and S4-3 



further confirm its thermo-mechanism. But, the crystal 
structure of r27RCL has not been reported before. Even 
if it is available, it is still difficult to predict the specific 
thermostabilizing mutations based on its structure be- 
cause there are so many structural features that stabilize 
proteins [37] and the contribution of each feature to the 
total free energy of thermostability is quite small. 

Conclusion 

P. pastoris proved to be a valuable cell factory for the 
directed evolution of the lipase from R. chinensis for 
thermostability improvement, especially when combined 
with the novel method of library construction described 



Enzyme 


K m (mmol/min/mg) 


'tea, (S 1 ) 


U« m CM " 1 s 1 ) 
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6.22 X10 4 


S4-3 
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here. In the present study, the thermostability of r27RCL 
was improved significantly by two rounds of error-prone 
PCR and two rounds of DNA shuffling in P. pastoris. 
The S4-3 variant is the most thermostable lipase from 
Rhizopus sp., under the conditions tested. Compared 
with the parent r27RCL, the optimum temperature of 
S4-3 was two degrees higher, T m was 22 degrees higher 
and half-lives at 60°C and 65°C were 46- and 23- times 
longer, respectively. Moreover, the catalytic efficiency 
k cst IK m of S4-3 was comparable to the parent. In sum- 
mary, S4-3 maintained good catalytic properties and 
exhibited much greater tolerability and stability at high 
temperatures, which made this variant a promising can- 
didate for industrial applications at high temperatures. 

Methods 

Strains, plasmid, reagents and cultures 

Restriction endonucleases, Taq DNA polymerase were 
purchased from Takara. All chemicals were of analytical 
grade or higher quality and purchased from Sigma. P. 
pastoris GS115 {His' Mut + ) and the pPIC9K expression 
vector are from Invitrogen BV. Recombinant plasmid 
pPlC9K-proRCL and strain GS115/ pPlC9K-proRCL 
were constructed by Yu et al. [10]. The media used in 
this study include YPD (lOg/L yeast extract, 20g/L pep- 
tone, 20g/L dextrose), MD (13.4g/L YNB, 0.016 uM 
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Figure 5 Thermostability of the parent and S4-3 at 60°C (a) and 65°C (b) 
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Figure 6 Three-dimensional structure of r27RCL. Locations of six 
amino acid substitutions (A129S, A230T, S234F, L180H, T218S, H317P) 
in S4-3 were shown in red sticks and the catalytic residues were 
shown in blue sticks (in the oval) 



biotin, 20g/L dextrose), BMGY (lOg/L yeast extract, 
20g/L peptone, lOg/L glycerol, pH 6.0 100 mM potas- 
sium phosphate, 13.4g/L YNB, 0.016 uM biotin), BMMY 
(lOg/L yeast extract, 20g/L peptone, lOg/L methanol, 
pH 6.0 100 mM potassium phosphate, 13.4g/L YNB, 
0.016uM biotin), BMMYA (lOg/L yeast extract, 20g/L 
peptone, 5g/L methanol, pH 6.0 100 mM potassium 
phosphate, 13.4g/L YNB, 0.016 uM biotin, 6g/L agar). 

Production of proRCL variants by error-prone PCR 

The first two rounds of random mutagenesis of proRCL 
gene were carried out by error-prone PCR as described 
elsewhere [43]. A 50 uL reaction mixture contained 5 u_L 
PCR Gold buffer, 0.7 mM MnCl 2 , 3 mM MgCl 2 , 0.5 mM 
dATP/dGTP, 2.5 mM dCTP/dTTP, 40 pmol of each pri- 
mer (BC-F and BC-R), 5 ng of plasmid pPlC9K-proRCL, 
and 2.5 U of Taq polymerase. PCR cycle conditions 
were: initial denaturation for 3 min at 94°C; followed by 
30 cycles of 30 s at 94°C, 30 s of annealing at 55°C and 
1 min of extension at 72°C; and a 10 min final extension 
at 72°C. The template for the second round of error- 
prone PCR used an equimolar mixture of positive var- 
iants from the first generation of random mutagenesis. 

Production of proRCL variants by DNA shuffling 

The last two rounds of random mutagenesis of proRCL 
gene were carried out by DNA shuffling based on the 
previously methodology reported by Zhao and Arnold 
[44]. The DNA template was an equimolar mixture of 
positive variants from the former round. The template 
was digested with 0.1 U DNase I in 200 uL of 50 mM 
Tris-HCl pH 7.4, 1 mM MgCl 2 for 10 min at 37°C. 



Fragments of 50-100 bp were purified by gel extraction 
(QIAGEN). Then purified fragments were resuspended 
in PCR mixture (without primers) containing 5 u_L PCR 
Gold buffer, 200 uM of each dNTP and 1.25 U of Taq 
polymerase. PCR cycle conditions were: initial denatur- 
ation for 3 min at 94°C; followed by 10 cycles of 30 s at 
94°C, 30 s of annealing at 45°C and 1 min of extension 
at 72°C; 10 cycles with annealing temperatures of 48°C 
and 51°C, respectively; and a 10 min final extension at 
72°C. Then, 1-3 u± of primerless PCR product was 
added to a 50 uL PCR mixture with 2 uM primers (BC- 
F/BC-R). Other PCR conditions were similar to those 
described above, and a single product comprising full- 
length proRCL was obtained. 

Production of his4 3'-end containing fragment (vector 
fragment A and D) 

Expression vector fragments were prepared by PCR 
amplification. The plasmid pPlC9K-proRCL was first 
amplified using flanking primers A-F/A-R or D-F/D-R to 
produce vector fragment A or D (Figure 1). The PCR 
was carried out in a final volume of 50 uL containing a 
mixture of 50 ng plasmid pPIC9K-/>ra/?CL, 1 uM of each 
primers, 5 uL PrimeSTAR DNA polymerase buffer, 
100 uM of each dNTP and 1.25 U of PrimeSTAR DNA 
polymerase. PCR cycle conditions were: initial denatur- 
ation for 3 min at 94°C; followed by 30 cycles of 30 s at 
94°C, 30 s of annealing at 50°C and 2.5 min or 7.5 min 
extension at 72°C for fragment D or A, respectively; and 
a 10 min final extension at 72°C. 

Library construction 

P. pastoris was transformed with 2 ug of fragments A, D 
and B/C in a 10:1:1 molar ratio by electroporation 
according to the protocol reported by Wu and Letch- 
worth [45]. Briefly, the competent cells, pretreated with 
0.1 M lithium acetate and 10 mM dithiothreitol, were 
transformed with 2 |ig DNA fragments by electropor- 
ation with instrument settings of 1.5 kV, 25 uF and 186 
€1. Transformed yeast colonies appeared in 4-6 days at 
30°C on MD plates. The variants were confirmed by 
DNA sequencing (Sangon Biotech Company). 

Overlapping sequence length 

In order to optimise the recombination efficiency, the 
length of the overlapping sequences between vector and 
gene fragments was optimized. It was noteworthy that 
the lipase gene fragment was be produced by normal 
PCR (called L) instead of ep-PCR. We generated frag- 
ment L using 10 bp, 25 bp, 40 bp, 60 bp, 70 bp and 
85 bp of overlapping primers listed in Table 4. These 
fragments were mixed with fragment A and D in a 
10:1:1 molar ratio and electroporated into P. pastoris. 
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Table 4 Oligonucleotides utilized in this study 



Primer Oligonucleotide sequence 

BC-F GCTGCrAAAGAAGAAGGGGTATCTCTCGAGAAAAGAGAGGCTGAAGCTTACGTAGAATTCCCTAGG 

BG-R GTAAGTGGCCAAGTTGAACTGAGGAACAGTCATGTCTAAGGGGAATTAA7TCGCGGGCGC 

A-F GGCCGCGAATTAATTCGCCTTAGACATG 

A-R AGATCTTGATATAAATrrCAGGTTTAAAATC 

D-F TCGACAATTGGTTTGACTAATrCCATAATCTG 

D-R ACGTTTCGTCTTTGGATGTTAGTCT 

BC-F 10 ATTCCCTAGG 

BC-R10 TCGCGGCCGC 

BC-F25 TGAAGCTTACGTAGAATTCCCTAGG 

BC-R25 CTAAGGCGAATTAATFCGCGGCCGC 

BC-F40 CGAGAAAAGAGAGGCTGAAGCTTACGTAGAATTCCCTAGG 

BC-R40 GAGGAACAGTCATGTCTAAGGCGAATTAATTCGCGGCCGC 

BC-F55 GAAGGGGTATCTCTCGAGAAAAGAGAGGCTGAAGCTTACGTAGAATTCCCTAGG 

BC-R55 GCCCAACTrGAACTGAGGAACAGTCATGTCTAAGGCGAATTAATFCGCGGCCGC 

BC-F70 CATTGCTGCTAAAGAAGAAGGGGTATCTCTCGAGAAAAGAGAGGCTGAAGCTTACGTAGAATTCCCTAGG 

BC-R70 CGGTCTTCTCGTAAGTGCCCAACTTGAACTGAGGAACAGTCATGTCTAAGGCGAATTAATTCGCGGCCGC 

BC-F85 TACTACTATTGCCAGCATTGCTGCTAAAGAAGAAGGGGTATCTCTCGAGAAAAGAGAGGCTGAAGCTTACGTAGAATTCCCTAGG 

BC-R85 AGAATCTAGCAAGACCGGTCTTCTCGTAAGTGCCCAACTTGAACTGAGGAACAGTCATGTCTAAGGCGAATTAATTCGCGGCCGC 

BC-F65 CTGCTAAAGAAGAAGGGGTATCTCTCGAGAAAAGAGAGGCTGAAGCTFACGTAGAATTCCCTAGG 

BC-R65 TTCTCGTAAGTGCCCAACrTGAACTGAGGAACAGTCATGTCTAAGGCGAATFAATFCGCGGCCGC 



The number in the name of the primers indicated the length of the overlapping sequences. 



Pre-screening 

Colonies from libraries of each generation were repli- 
cated with sterile toothpicks from MD plates onto 
YPD plates (for the storage of variants) and BMMYA 
plates (for gene expression and screening), simultan- 
eously. The parent strain was used as a positive con- 
trol. YPD and BMMYA plates were then incubated at 
28-30°C for 2 d and 5 d, respectively. Then, BMMYA 
plates from each ep-PCR and DNA shuffling gener- 
ation were subjected to heat-treatment at 45°C, 50°C, 
55°C and 65°C, respectively, for 60 min. Screening was 
performed with 15 mL of top agar (5g/L) containing 
160 uL of Fast blue RR (89g/L in dimethyl sulfoxide) 
and 80 uL a-naphthyl acetate (40g/L in dimethyl for- 
mamide) [46]. Positive colonies exhibited dark brown 
coloration within two minutes after pouring the top 
agar and were further screened using the 96 MTPs 
screening method. 

Screening for evolved variants 

All positive candidates screened from top-agar staining 
were picked and cultivated in deep-well microtiter plates 
containing 900 uL of BMGY medium, at 28°C, 300 rpm 
for 16 h to reach an OD 600 of 1. Then, cultures were 
harvested by centrifugation at 3000 g at 4°C for 10 min 



and the cells were suspended in 300 uL of BMMY 
medium for gene expression. The cultures were kept at 
28°C, 300 rpm for 84 h with supply of 100 uL fresh 
BMMY medium and lOg/L methanol every 24 h. 

Each culture was centrifuged at 5000g at 4°C for 
10 min. The supernatant was then transferred into 96- 
well microtiter plates for heat-treatment at 60°C or 65°C 
for 60 min. Initial and residual lipase activities were 
measured. All the assays were done in triplicate and sig- 
nificant differences (p < 0.05) were measured. Only var- 
iants that showed both higher ratio of residual activity 
and similar initial activity compared to the parent en- 
zyme were selected. 

Enzyme purification 

A single colony from each variant was cultivated in 
100 mL of BMGY medium shaken at 28°C and 250 rpm 
in 500 mL glass flasks. When cultures reached an OD 600 
of 1, the cells were centrifuged and resuspended in 
25 mL of BMMY medium to obtain an QD 60 o of 4 and 
shaken at 28°C and 250 rpm in 250 mL glass flasks for 
84 h. The cultures were supplemented with methanol 
(5g/L) to induce the expression of lipase every 12 h. The 
culture was then centrifuged and the supernatant was 
collected for protein purification. 



Yu et al. Microbial Cell Factories 2012, 11:102 
http://www.microbialcellfactories.eom/content/1 1 /1 /1 02 



Page 1 0 of 1 1 



Cell free medium was concentrated and interchanged 
with 10 mM Tris-HCl buffer (pH 7.5) by ultra filtration 
through a 10-kDa membrane (Millipore). The concen- 
trated solution was loaded onto a SP-Sepharose column 
(AKTA, 1 cm x 20 cm) equilibrated with 20 mM Tris- 
HCl buffer (pH 7.5) and eluted with 0-0.5 M NaCl in 
the same buffer. Fractions containing lipase activity were 
pooled, concentrated and loaded on a Phenyl-sepharose 
6 FF column (AKTA, 1.6 cm x 20 cm) equilibrated in 
50 mM Tris-HCl buffer (pH 7.5) containing 1.6 M am- 
monium sulfate. Lipase was then eluted in an ammo- 
nium sulfate concentration gradient decreasing from 1.6 
to 0 M in 50 mM Tris-HCl buffer (pH 7.5) and 4 mL 
fractions were collected at a flow rate of 0.8 mL/min. 

Lipase properties 

Lipase activity was measured on emulsified />-nitrophe- 
nyl palmitate (pNPP) according to Kordel et al. [47]. 
Unless stated otherwise, the typical enzymatic reaction 
was carried out at 40°C, pH 8.5. One enzyme unit was 
defined as the amount of enzyme releasing 1 umol of p- 
nitrophenol per minute under the assay conditions. Op- 
timal temperature (T opt ) was determined by measuring 
the enzyme activity at pH 8.5 under various tempera- 
tures (20-60°C). T m values were determined as 
described previously [43]. 10 U purified enzyme solu- 
tions were incubated for 30 min in the temperature 
range of 20-80°C, cooled in the ice bath for 20 min and 
equilibrated at room temperature for 5 min. Then, re- 
sidual activity was determined at 40°C and expressed as 
a percentage of the initial activity. The temperature at 
which 50% of lipase activity lost was the T m value. Puri- 
fied enzymes were incubated at specific temperatures 
(60°C and 65°C) for half-life (ti/ 2 ) determination. 10 U 
purified enzyme solutions were taken at various time 
intervals, cooled in the ice bath for 20 min and equili- 
brated at room temperature for 5 min. Then, residual 
activity was determined and the incubation time at 
which 50% of lipase activity lost was the t 1/2 value. In- 
activation was followed until > 80% of the activity was 
lost. The Michaelis-Menten kinetic parameters /c cat 
and K m were calculated using £>NPP as substrate. 
Lineweaver-Burk plots were used to determine k cat and 
K m parameters, assuming that the reactions followed a 
simple Michaelis-Menten kinetics. 

Thermo-mechanism analysis 

A three-dimensional model of R. chinensis lipase was 
built by SWISS -MODEL protein automated modelling 
program [48] on the basis of crystal structure of 1LGY 
(crystal structure of lipase II from Rhizopus niveus 
solved with a resolution of 2.20 A) [49], which showed 
the highest homology of 80.38% to R. chinensis lipase. 
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